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Abstract. We review the properties of biological motor proteins which move along 
linear filaments that are polar and periodic. The physics of the operation of such motors 
can be described by simple stochastic models which are coupled to a chemical reaction. 
We analyze the essential features of force and motion generation and discuss the general 
properties of single motors in the framework of two-state models. Systems which contain 
large numbers of motors such as muscles and flagella motivate the study of many 
interacting motors within the framework of simple models. In this case, collective effects 
can lead to new types of behaviors such as dynamic instabilities of the steady states 
and oscillatory motion. 



1 Introduction 

An essential and striking features of living cells is their ability to generate mo- 
tion and forces. Important examples are cell motility, the contraction of muscles 
but also active phenomena within cells that allow for material transport and the 
motion of organelles, for example during cell division and mitose. These move- 
ments and forces are generated on the molecular level by protein molecules that 
are driven by chemical reactions in a far from equilibrium situation. Important 
examples are motor proteins, enzymes which are specialized to work as motors. 
In eucariotic cells, several families of motor proteins exist which interact with the 
cytoskeleton, a complex three-dimensional elastic network of long rod-like fila- 
ments inside the cell which is essential for the mechanical stability and integrity 
of the cell [1]. 

A motor protein of the cytoskeleton interacts specifically with a certain type 
of filament along which it is able to move in presence of Adenosinetriphosphate 
(ATP) which is a chemical fuel. The filaments serve as guides or tracks for the 
motion. Two types of filaments play this role: microtubules and actin filaments. 
Both are formed by a polymerization process from identical monomers (actin 
and tubulin monomers, respectively), leading to a regular and periodic struc- 
ture. An important feature is their polarity: The filaments are asymmetric with 
respect to their two ends. This symmetry has its origin in the asymmetry of the 
monomers which form a polar filament structure with two different ends which 
are denoted "plus end" and "minus end" . This polar symmetry is essential for 
motor operation as it defines the direction of motion. Motor proteins are clas- 
sified into several families: myosins, kinesins and dyneins. Myosins move always 
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along actin filaments and towards the plus end. Kinesins and dyneins move along 
microtubules, kinesins move towards the plus end and dyneins towards the minus 
end, see Fig. 1. 




Fig. 1. Schematic representation of molecular motors and track filaments, (a) myosin 
interacting with actin filaments, (b) kinesin moving along a microtubule. Both types 
of filaments are polar and periodic, their two different ends are denoted "plus" and 
"minus" . 



Myosins are prominent for their role in the contraction of muscles [2, 3]. In 
this case, many myosin molecules form a linear structure, a myosin filament, 
which interacts with actin filaments arranged in parallel. The action of myosin 
motors then induces the relative sliding of the two types of filaments. In muscle 
cells, a very large number of filaments is organized in a regular way which on a 
macroscopic scale leads to the muscular contraction. Myosins also occur within 
normal cells where they play an important role for cell motility and the organi- 
zation of actin. Kinesins occur in large numbers in neurons, where they play a 
key role in transport of vesicles along the axon towards the synapse. Both types 
of motors have two identical heads of a size of about 10-20nm which are the 
elementary force-generating elements as well as a tail which is used to attach 
the motor to another structure [4]. 

The energy source of this process is the hydrolysis reaction ATP ADP+P 
of ATP to ADP and Phosphate (P). The motor protein M (or more precisely, 
the head domain containing the ATP binding site) undergoes a chemical cycle: it 
binds ATP and hydrolizes the bound ATP. Subsequently it releases the products 
ADP and P. We denote the different chemical states by M, M-ATP, M-ADP-P, 
and M-ADP, respectively. After completion of the cycle (M-I-ATP ^ M-ATP 
M-ADP-P^ M-ADP+P M+ADP+P) the motor is unchanged. However, dur- 
ing this process it has undergone conformational changes and it has hydrolyzed 
one ATP molecule, see Fig. 2. The different conformations which occur during 
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Fig. 2. Chemical cycle of a motor molecule M. After completion of the cycle on 
molecule ATP is hydrolyzed to ADP and phosphate (P). 



the chemical cycle in general have different geometries and properties and can in 
particular have different interaction characteristics with respect to the filament 
[5]. As a result, the motor protein undergoes chemistry-driven changes between 
strongly and more weakly bound states ("attachments" and "detachments"). 
This coupling between chemistry and binding permits the creation of motion 
along a polar filament [6, 7, 8]. 

The force and motion generation of individual or groups of motor proteins 
can be studied experimentally by a variety of techniques. In so-called motility 
assays, motors are attached to a substrate [9, 10, 11, 12]. Filaments in solution 
bind to the motors and in presence of ATP start moving along the surface. 
In these experiments typically several or many motors interact with a single 
filament. 

In order to observe the forces generated by an individual motor, the pro- 
cessivity of the motor becomes important. Myosin is not processive. During the 
chemical cycle it detaches from the filament during a significant period of time. 
During this time it can easily diffuse away from the filament if it is not held 
in place. Forces generated by single myosin molecules have been observcid by 
different techniques, see Fig. 3 (a) and (b). Micron sized bead have been coated 
with low density of motors and optical traps have been used to bring a filament 
in contact with the bead and possibly only a single motor [13]. Forces on the 
filament are then measured by observing the displacement in the trap. Another 
possibility is to fix the actin and to manipulate the bead with an optical trap. 
Such experiments reveal that the motor induces stochastic displacements of the 
order of 5 — lOnm which last for several milliseconds and peak forces of the order 
of IpN [11, 13]. 

Kinesin, is a processive motor [14]. A single kinesin moves along microtubules 
for a distance of the order of micrometers before losing the filament. Motion of 
single kinesin molecules can therefore be observed directly, either by attaching 



4 Frank Jiilicher 




myosin 



Fig. 3. Examples for micro-manipulation experiments, (a) Forces induced my myosin 
on actin can be measured by the deflection of a micro-needle [If], (b) Force measure- 
ment using optical traps [13] . (c) Forces generated by a single kinesin molecule observed 
by displacements of a bead in an optical trap [15]. (d) Forces induced by an electric 
field i5 in a motility assay using linear grooves to orient the filaments [18]. 



a small bead to the motor and observing the displacement of the bead or by 
directly marking the molecule with a fluorescent dye. Optical traps have been 
used to study the velocity of motion as a function of an appHed load [15]. It 
has been shown, that kinesin moves in a step-wise fashion with characteristic 
steps of 8nm size [16]. This step size coincides with the period of kinesin binding 
sites along microtubules which demonstrates that the ATP-driven reaction cycle 
induces steps between periodically spaced binding sites. The characteristic time- 
scales for kinesin motion are milliseconds and forces up to 5pN are generated 
[15, 17]. 

The standard explanation for the processivity of Kinesin is based on the 
fact that each molecule consists of two head domains which both hydrolize ATP 
and undergo the chemical cycle in a coordinated way [14, 19]. In this situation, 
processivity is possible if both heads very rarely detach from the filament at 
the same time thus allowing the motor to keep attached while displacing. Single 
heads do have the capability to generate forces and motion but in most cases are 
not processive. Recent experiments however show evidence that single kinesin 
heads can move processively in certain cases [20]. Myosin molecules also often 
have two heads. However, even with two heads they do not become processive 
since every myosin head has a tendency under normal operating conditions to 
be unbound during about 90% of its chemical cycle. 

From the point of view of a general classification of the physical mechanisms 
which can lead to motion generation, motor proteins fall in a class of systems 
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characterized by the fact that they operate on molecular scales and generate 
motion along a one- dimensional polar and periodic structure (sometimes called 
ratchet [21]) by a non-equilibrium rectification process [22, 23, 24, 25, 26, 27, 
28, 29, 30]. For a system on a molecular scale, fluctuations play an important 
role for its function and properties. These fluctuations can be both of thermal 
origin or they can arise due to the stochasticity of individual molecule chemical 
reactions. Therefore, a theoretical description requires the use of concepts of non- 
equilibrium statistical physics. Biological motor proteins are the most prominent 
examples for these systems. However, there is a growing number of artificially 
designed "molecular motors" [31, 32, 33, 34], suggesting that the physics of these 
systems is relevant for micro- and nano-technological devices. 

In the subsequent sections, we describe a simple modelization of the generic 
aspects of these systems, using biological motor proteins as a guide. In section 
2 we introduce the basic concepts of such a modelization. Properties of single 
motor motion which follow from this description are described in Sect. 3. In 
section 4, we discuss the conscqiiences of many motor systems and demonstrate 
that new phenomena such as dynamical instabilities follow naturally from these 
models. Finally, we present in section 5 a discussion and an outlook. 

2 Simplified models 

In order to keep the description simple and to focus on generic properties, we do 
not aim to capture microscopic structural details of biological motor proteins. 

We use a simplified picture taking into account the periodicity and polarity of the 
linear track and the fact that during the chemical cycle conformational changes 
occur [8]. 

The main assumption required for this simplification is a separation of time- 
scales [6, 8]. Even though a macromolecule like a protein is characterized by 
a large number of microscopic degrees of freedom, most of them relax on time 
scales shorter than the typical relevant time scales for the chemical cycle. These 
degrees of freedom are therefore to a good approximation thermally equilibrated. 
Only a few slow degrees of freedom have to be described by dynamical equations. 
These relevant degrees of freedom arc collective modes of the system. Examples 
are chemical reaction coordinates but also the overall position variable of the 
motor with respect to the filament. The latter is important in order to describe 
the coupling between chemistry and motion. 

2.1 Energy landscapes and chemical transition rates 

Following Ref. [6] , we assume for simplicity that the chemistry of a single head 
can be described by a number m of discrete states or conformations i = l..m. 
Many biochemical models focus on the four states M, M-ATP, M-ADP-P, M- 
ADP but other reaction intermediates could also be included [6, 35]. Transitions 
between these states are fast compared to the typical time between transitions 
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and the mechanical action or motion of the motor. Therefore we use a chemical 
kinetics description for changes between states. 

If we consider a motor in conformation i, we can define a potential or inter- 
action energy profile along the filament. Suppose that one small region of the 
motor, e.g. in the tail, is used to transmit forces or to attach a cargo. We imag- 
ine this point to be held at a position x along the filament. We can now define 
Wi{x) to be the energy of the motor, including possibly bound ATP, ADP or P, 
and including the energy of the filament as the motor is held at position x. This 
total energy is in fact an effective free energy defined formally by integrating 
over all rapidly relaxing microscopic degrees of freedom but keeping the enzyme 
in its chemical state. The conformation of the system motor-filament is then 
fully characterized by the pair {i,x} of internal state and position with respect 
to the filament. Note, that the potentials refiect the symmetry properties of the 
filament. If the filament is polar and a periodic array of identical monomers, 
the potentials are periodic with period I, Wi{x) = Wi{x + I) and asymmetric, 
Wi{x) ^ W,i-x). 

In order to describe the dynamics of the system {i, x}{t), we use a stochastic 
overdamped dynamics at constant temperature T within a given state i 

rji^^x = -d,Wi{x) + Q{t) . (1) 

Here rji is a protein friction and ((t) is a Gaussian white noise in state i with 
< Ci{t)Cj{0) >= 2?7i^(t). The chemical reactions between states {i,x} and {j,x} 

U)ij{x) 

{i,x} # {j,x} (2) 

UJji{x) 

occur with Poisson statistics with reaction rates LUij{x). Since the position vari- 
able X is also a conformational degree of freedom (the motor in general changes 
its shape while displacing), transition rates are in general x-dependent. Note, 
that for simplicity in Eq. (2) we have assumed that transitions between states 
happen instantaneously and without displacement. Furthermore, we have used 
the fact that thermal relaxation is very fast compared to the chemical cycle and 
all rapid degrees of freedom are equilibrated at constant temperature T. In fact, 
the typical relaxation time of temperature gradients which have developed on a 
length scale I can be estimated as r = CI'^/k, where C is the specific heat of 
the material per volume and k the thermal conductivity. Using typical values for 
water and length scales of the order of lOnm we find r ~ 10~^ — 10~*s, which 
is fast compared to typical cycle times of several ms. This argument shows that 
the motor operates isothermally, i.e. temperature gradients are not created and 
cannot be used to generate motion as e.g. in the case of Feynman's ratchet [21]. 

It is now convenient to use a Fokker-Planck description [36] and to introduce 
distribution functions Pi{x,t) for the probability to find within an ensemble of 
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identical systems the motor at time t at position x in state i. These distributions 
then obey the equations 

dfPi + dxJi = '^{t^3i{x)Pi{x) - ujij{x)Pj{x)) (3) 
Ji = r,i\-kBTd^Pi-PidxWi + PiU^) . (4) 
The total density and total current 

P{x,t) = Y.Pi (5) 

1=1 

m 

J{x,t) = Y,Ji (6) 
j=i 

obey the conservation law dtP+dxJ = 0. The average velocity in the steady state 
with stationary and periodic distribution function Pi{x) = Pi{x + l), dtPi = is 
given by 

v = I Jdx/ I Pdx (7) 
Jo Jo 

In order to characterize the chemical rates, we first introduce the chemical 
potentials of the fuel and hydrolysis products in bulk solution. We denote hatp^ 
fJbADP and jjLp the free energy per ATP, ADP or P molecule, respectively. As 
an illustrative example, we first consider the four chemical states M [i = 1), 
M-ATP (i = 2), M-ADP-P {i = 3) and M-ADP (i = 4), often encountered for 
biological motor proteins. A general reaction kinetics for all eight reaction rates 
which is consistent with the ATP hydrolysis reaction can be written by using 
four different kinetic coefficients: 

wi2 = ai exp[(Wi + iiATp)/kBT] W21 = ai eyoi,\W2/kBT] 

W23 = Oi2 exp[M^2/fcBT] W32 = Q!2 exp[W3 / ksT] 

W34 = as expliys/fcsT] W43 = 03 expi(W4 + ^ip)/kBT] ' 

W41 = Q;4exp[(W4 + Hp)/kBT] LUii = a4exp[(T4^i + ^adp + tJ'p)/kBT] 

(8) 

Here, we have used the condition of detailed balance of the rates. The functions 
ai{x) characterize the possible reaction scenarios. Note, that since transitions 
are fast and therefore occur for fixed x, the c;licniical rates do not depend on 
the external force /ext or local stresses. Force-dependent are only displacements 
described by equation (1). The present modelization differs in this respect from 
models which use discrete transitions also to describe displacements. In this case 
chemical rates are strain dependent [6, 37, 43]. From a physical point of view, 
all models are of course equivalent. 



8 Frank Jiilicher 



2.2 Two-state models 

It is useful to further simplify the generic description introduced above. The 
m-state model allows in principle to describe many details of the chemical cycle 
and the various conformations of the motor. However, it contains a large number 
of free parameters which are unknown. It is therefore useful to further simplify 
the model. In fact, in order to describe physical aspects of motion generation 
and force generation, it is sufficient to keep only two different states [24, 26, 27]. 




X 

Fig. 4. Two state model defined by two polar and periodic potentials Wi and W2 as 
well as periodic transition rates coi and lj2- Pumping between the two states induces 
average motion. 

We rewrite the Fokker-Planck equations (4) for two states i — 1,2: 

dtPl + d^Ji = -UJi{x)Pi + U02{X)P2 

dtP2 + d.,J2 = LOi{x)Pi - W2{X)P2 , (9) 

where we have introduced lji ~ UJ12 and uj2 = LO21 and the currents are the same 
as introduced in Eq. (4). This system is sketched in Fig. 4 for an example of 
shifted periodic and asymmetric potentials. 

This two-state model is still very flexible and allows to describe situations 
which capture many of the physical aspects of biological protein motors. Fig. 5 
shows choice of potentials Wi and W2 adapted to the commonly accepted pic- 
ture of myosin function [7] : a myosin head detaches from the actin filament after 
binding ATP. In the unbound state ATP is hydrolyzed (M-ATP^ M-ADP-P). 
The head (M-ADP-P) is now again able to bind actin. As it encounters a binding 
site along the filament, it re-attaches under phosphate release. After reattach- 
ment, a force-generating step occurs and ADP is released, which completes the 
chemical cycle. As illustrated in Fig. 5 this process can be captured by two dif- 
ferent potentials, Wi and W2 representing the unbound state (a flat potential) 
and the bound state (a potential with periodic structure) respectively. After 
the force-generating step (power-stroke), the x variable has reached a potential 
minimum. Here, the system can be actively excited to the unbound state under 
ATP binding. As it reattaches to the filament, the slope of the potential refiects 
the mechanical force generated at this point. A displacement is generated as the 
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Fig. 5. Two state model representing a situation motivated by the functioning of 
myosin. (1) Binding an ATP molecule, myosin detaches from actin (2). After hydroly- 
sis (3), it rebinds and generates a force (4) and a displacement. In a two state model, 
two potentials Wi and W2 characterize attached and detached states with the tail at 
position X. The shaded areas are "active sites" where ATP-driven transitions occur. 



system slides downhill along the energy profile to reach the potential minimum. 
Note, that microscopically this displacement could correspond either to a tilt of 
the head domain as sketched in the figure or to other more complex processes. 
The microscopic structure associated with this displacement is not characterized 
by this description. 

In a two-state picture, the chemical reaction cycle as described by the kinetic 
equations (8) has to be divided into two substeps. One possibility is to introduce 
the forward and backward rates ai and a2 for the combined process of ATP- 
binding and hydrolysis 

ai 

M + ATP ^ M-ADP-P , (10) 
a2 

and the rates /3i and P2 which describe the process of product release and bind- 
ing: 

Pi 

M + ADP + P ^ M-ADP-P . (11) 

P2 
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The complete chemical cycle is now the subsequent transitions ai and (32- As 
long as a2 and j3\ are nonzero, there is a nonvanishing probability for an inversion 
of the cycle (i.e. ATP generation) by following the steps and fii. We define 
Wi to be the energy of a free motor together with the product molecules (M + 
ADP + P) and to be the energy of the complex M-ADP-P. Detailed balance 
of the chemical reactions then implies 



"1 ^ ^{Wi-W2+A,j.)/kBT (^-^2) 
0.2 



/32 

where we have introduced the chemical driving force 

All = ^lATP - MADP - MP • (14) 

The transition rates of the two-state model are the superpositions = aj + /Jj. 
Introducing two unknown functions a{x) and (5{x) which describe conformation 
dependent energy barriers, we can therefore write 

uj^{x) = a(a;)e(^^+^'')/'=-^ + /3(a;)e^^/'=-^ 

u;2{x) = [a{x)+(3{x)]e^^/'"'^ . (15) 

Note, that other choices to divide the reaction cycle in two relevant substeps 
leads to the same result, but redefines the arbitrary functions a and (3 and shifts 
the potential W2 by a constant value. 

The functions a{x) and (3{x) define the coupling of the chemical reaction 
to conformation. Very important is the concept of localized or conformation- 
dependent transitions where the functions are peaked within a narrow a;-interval 
but negligible outside this interval. An example is the ATP-binding step which 
in Fig. 5 is restricted to occur within an "active region" of conformation space 
corresponding to the potential minimum while the conformations at the begin- 
ning of a force-generating power stroke are not supposed to bind ATP. As we 
will describe in the subsequent sections, the localization of transitions via the 
functions a and /? plays an important role for many interesting cases. 

Similar to the case of myosin, the two state model can also be adapted to 
other situations such as the motion of kinesin molecules containing two heads 
both of which contain an ATP binding site. In principle, a general description 
would require eight internal states and a complex reaction scenario. A possible 
simplification arises from the idea of a coordinated binding and unbinding of 
the two heads in a hand-over-hand fashion as shown schematically in Fig. 6 [14]. 
Such a coupling would reduce the number of relevant degrees of freedom [38]. In a 
two-state model, this feature can be captured in a simplified way by associating 
each state with one of the heads being bound to the filament. Denoting the 
motor with head 1 or head 2 bound to the filament by Mi and M2, respectively, 
we define the energies of these situations Wi and W2. Because the two heads 
forming a kinesin motor are identical, the energy landscapes of both states are 



Force and Motion Generation of Molecular Motors: A Generic Description 



11 








2 




7 




X 



Fig. 6. Hand over hand motion suggested for kinesin. At any time, one head is bound 
and the second head moves towards the next binding site. This situation can be repre- 
sented by describing both heads by identical but shifted energy landscapes Wi (x) and 

W2{x). 



two identical potential profiles which are shifted with respect to each other by one 
monomer period 1/2 on the filament: Wi{x) = Wi{x + l); W2{x) = Wi{x + 1/2), 
see Fig. 6. However, the characteristic step-size of an individual head I in this 
picture corresponds to two monomer sizes. Therefore, the potentials Wi are l- 
periodic, while the total system is invariant under x —^ a: -I- //2 if at the same time 
the two states are exchanged. Therefore the transition rates also obey lui{x) = 
uji{x + l) and 0^2(2;) — uJi{x + l/2). Again, the idea of active regions and localized 
transitions is important. Assuming that all transitions occur at conformations 
which correspond to the potential minimum, we obtain a system which operates 
in an almost deterministic way where the chemical cycle is closely correlated to 
a particular displacement. In the case of non-localized transitions, the chemical 
cycle is related to motion in a more irregular way. 

A useful representation of the transition rates in the hand-over-hand picture 
is to assume that a full ATP hydrolysis cycle changes Mi to AI2: 

Ml + ATP ^ M2 + ADP + P . (16) 

Because of the symmetry between the two heads, the reaction M2 + ATP ^ 
Ml + ADP + P occurs with the same rates. This leads to the total transition 
rates 

LUiix) = e'^i(-)/fcBT[^(^)g4M/fcBT ^ ^ 1^2)] 
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The unknown function a{x) = a{x + 1/2) is /-periodic. Note, that the clioice 
given in_Eq. (17) is a special case of Eq. (15) if we identity Afj. = 2Aij,, a{x) = 
a(x)e-^''/'=BT^ ^ ^(2; + l),W2 = W2+Aii and Wi = Wi. This example 
demonstrates that Eq. (15) is a general choice which can describe very different 
types of couplings of an ATP hydrolysis cycle to a two state model. 

We have now defined a two-state model which can describe a system un- 
dergoing an ATP hydrolysis cycle and moving along a periodic structure. As 
discussed above, this model is sufficiently flexible to be adapted to situations 
which resemble the widely accepted pictures of myosin and kinesin functioning. 

3 Single motors 

We will now discuss general properties of the two-state model for a single motor 
introduced above [26, 39, 40]. Two generalized forces act on the system leading to 
an out-of equilibrium situation. These are the chemical "force" introduced 
in Eq. (14) and the mechanical force fext- If both generalized forces are kept 
constant, the system eventually attains a steady state with dfPi = 0. The steady 
state distribution functions satisfy two coupled differential equations of second 
order 



keTdlPi + {d^Pi){d^Wi - /ext) - PidlWi = r^icviPi - CJ2P2) 
kBTdlP2 + {d^P2){d^W2 - /cxt) - P2dlW2 = -vi^iPi - UJ2P2) , (18) 



where we have for simplicity assumes that the friction rj is the same for both 
states. This set of equations together with periodic boundary conditions Pi{x) = 
Pi{x + l) defines the steady state distributions. They can be calculated in special 
cases analytically, but in general numerical integration techniques are used. For 
each pair (zi/x, /ext), there is a uniquely defined average velocity 



Similarly, we can introduce the ATP hydrolysis rate r which denotes the number 
of chemical cycles performed per unit time [8]. Using the rates introduced in 
Eqns. (10) and (11), 




(19) 



where the satisfy the normalization condition 




(20) 




(21) 



If An = 0, the transition rates defined by Eq. (15) satisfy 



LU1/UJ2 = e 



-AW/ksT 



(22) 
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where 

AW{x) = W2{x) -Wi{x) . (23) 

This condition of detailed balance for the total transition rates indicates that 
transitions are just thermal fluctuations and that the system is not driven chemi- 
cally. If the external force also vanishes, the steady state is a thermal equilibrium 
with Pi = iVe~^'(^)/''^^ for which v = and r = 0. For Z\/i > 0, the system 
is chemically driven. If no external force is applied spontaneous motion with 
V ^ Q can occur, however only if the system is polar. For a symmetric system 
with Wi{x) = Wi{—x) and LLii{x) = uji{—x) the steady state distributions are 
also symmetric Pi{x) = Pi{—x). Since d^Wi is antisymmetric in this case, v = 
by symmetry according to Eq. (19). On the other hand, r is in general nonzero 
in this case (the functions ai are symmetric). For spontaneous motion to occur 
two requirements have to be fulfilled: detailed balance of the transition has to be 
broken, which corresponds to Zi/x ^ and the system must have polar symmetry. 
In the case of motor proteins the polar filaments play this role. 

In the presence of an external force /ext) the system can perform mechanical 
work, i.e. it operates as a motor. The work performed per unit time against the 
external force is 

W = -Utv (24) 
while the chemical energy consumed per unit time is given by 

Q = rAij, . (25) 

We can therefore define the efficiency of energy transduction 

This quantity is useful for forces applied opposite to the direction of motion 
where < 77 < 1. Note, that this definition relies on the fact that a bulk solution 
exists which plays the role of a thermodynamic reservoir and allows to define the 
chemical potential difference of fuel and products. In situations where reservoirs 
are small the efficiency would be more difficult to define. 

For a discussion of physical aspects of motion, it is useful to write 

uji{x)=uj2{x){^2{x)+e-'^'^"'-^) , (27) 

where 

Q{x) = ^-^^I^BT^^A^^Ik^T _ ^ (28) 

measures locally the rate of transitions violating detailed balance. Using the de- 
pendence of the chemical potential on particle concentration, /Xj = /x^-t-fc^T In Cj, 
we observe that 

^ / Catp ,n\ 
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where = e'^^'ATp-ADp-t^°p)/^BT ^j^g equilibrium constant of the hydrolysis 
reaction. Q therefore is a direct measure of the distance from chemical equilib- 
rium. Prom Eqns. (28) and (29) we find 

^^fzl/z for z\MABr«i 

\CATplCADpCpior Afx/kBT:$>l ^""^^ 

For our discussion of the two-state model it is useful to characterize the system 
by the functions L02 and i? instead of a and (3 which allows us to discuss motion 
generation without the need to introduce chemistry. This choice is more general 
and can be used also for cases where transitions between states are triggered 
by other processes than chemistry such as in artificially constructed systems 
[31, 32, 33]. 

The properties of this two-state models have been discussed in Refs. [8, 24, 
39, 40]. Calculating the average velocity w as a function of the externally applied 
force /ext typically leads to a behavior which is well approximated by a linear 
dependence 

V ~ ^;o(l - (/ext//«)) (31) 

for many different choices of the potential shapes and the transition rates. Here 
vo is the spontaneous velocity at zero force /ext = and the stalling force, 
i.e. the force for which the system stops moving. Deviations from this linear 
behavior mainly occur for forces larger than the stalling force |/ext| > |/s| or for 
forces parallel to its natural direction of motion fext/fs < 0. 

The observed force- velocity curves for kinesin motors show an almost linear 
behavior which can be characterized by vq and fs defined in (31). While fs ~ 
5pN does not depend much on experimental conditions, the no-load velocity vq 
depends on ATP concentration and attached viscous loads and is of the order of 
1/x/s or smaller [17, 41]. 

The orders of magnitude observed for kinesin can be reproduced by the two- 
state model. Using e.g. a choice of potentials as shown in Fig. 6 with transitions 
localized at the potential minimum, the stall force is approximatively given by 
the potential slope. Choosing a potential amplitude of ?7 ~ lOksT which is set 
by the available chemical energy of ~ 10 — IbksT and a period of I ~ 8nm 
of microtubules, this force is fs — U/l = 5pN consistent with the observed 
value. The spontaneous velocity of the two-state model can be estimated by 
Vo — l/{tc + tg), where tc is the time of the chemical steps and tg — Pr]/U 
is the sliding time in the potential. Therefore, the observation of vq does not 
fix both the chemical rate and the value of i] corresponding to protein friction. 
One estimate for the unknown friction rj is to assume a hydrodynamic friction 
with a viscosity rj^is a factor 10^ — 10^ larger than the one of water, suggesting 
V ~ ^vis' ~ (1 — 10)10~*kg/s. This value takes into account that the dissipation of 
proteins should be better represented by the viscous behavior of macromolecular 
solutions. A different approach is to assume that for large ATP concentration 
tc — uj2^il~^ is negligible and the friction 77 determines the sliding velocity 
Vo ~ U/lij. Estimating the maximal velocity to be Vynax ~ 10~^m/s, we find 
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T] ^ lO"'^ which can be seen as an upper bound since chemical steps which in 
general also contribute to friction are neglected. 

A key parameter characterizing the conditions of operation of the two-state 
model is the the dimcnsionlcss value U / QijJ2 which compares the typical chem- 
ical transition time with sliding times in the potential slope. With the arguments 
given above we estimate U/^Pf2LJ2 — 0.1 — 1, where we have used f2uj2 — 10^s~^ 
which is a typical transition rate [42]. However, different values arc also consis- 
tent with the observed force-velocity relation as the spontaneous velocity Vq is 
determined by the longest of the two time scales mentioned above. Additional 
information such as velocity fluctuations would be required to determine this 
value from experimental observations and to fix the orders of magnitude of all 
parameters of the model. 

The two state model is consistent with the observed behaviors for biologi- 
cal motor molecules and reproduces typical velocities and forces and the force 
velocity relation. Other types of models which use different representations of 
states and transitions have also be used to consistently describe the force- velocity 
relation of kinesin [38, 43, 44]. 

4 Collective Effects: Dynamic Instabilities 

In many biological situations, motor molecules and filaments do not operate as 
isolated enzymes but many motors are integrated in larger structures. Typical 
examples are actin/myosin in muscles and dyneins/microtubules in flagella and 
cilia. Furthermore, the presence on motor molecules in the cytoskeleton leads to 
complex physical properties of these systems on large scales [45] . 

The most prominent many-motor system is actin/myosin in muscles. In this 
case, myosin molecules are attached together by their tails to form a linear fila- 
ment. Myosin filaments and actin filaments are arranged in parallel in a highly 
organized fashion. In the presence of ATP they slide with respect to each other 
which macroscopically leads to muscle contraction. Experimental in vitro "motil- 
ity assays" can be used to study myosin function in an artificial environment. In 
these systems, myosin molecules are attached to a solid substrate using specific 
antibodies. Actin filaments in solution adsorb to the myosin coated surface and 
start to move in presence of ATP as a result of the action of myosin motors 
[7, 10, 12, 18]. 

Another example are cilia and flagella which are elastic linear extensions of 
many cells which generate a beating motion used to propel the cell within a 
solvent or to move the solvent. A flagellum typically contains 9 pairs of mi- 
crotubules, each pair coupled by a large number of dynein motors and other 
proteins which serve as structural elements. The motors create forces that lead 
to the bending and motion of the flagella. Interestingly, these motors are used 
to generate oscillating motion [1]. 

These systems demonstrate that the behavior of systems containing motors 
can on larger scales have new and different types of behaviors than the one ob- 
served for individual motor molecules. As a first approach to discuss the behavior 




of many motor systems, we generalize the two-state model to describe a large 
number of coupled motors [8, 46, 47]. A simple modelization of such a situation is 
sketched in Fig. 7 (a). Many motors which all are described by a two-state model 
are attached along a linear backbone with constant spacing s. Assuming that 
the spacing between motors is fixed implies that the backbone is rigid. Within 
this assumption all motors have the same velocity. In practical cases, elastic 
properties of the filaments and of a backbone coupling the motors can become 
important. For example in the case of muscles, the passive elastic behavior of 
proteins such as titin can play the role of elastically coupling motor-filament 
systems to their environment [1]. Similarly, in flagella, bending elasticity of the 
filaments is essential to allow for the generation of beating motion. Fig. 7 (b) 
and (c) sketches two simple ways to incorporate the effects of material elasticity 
in the modelization. 
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4.1 Mean field limit 

Motors coupled via a rigid backbone allow to illustrate the appearance of collec- 
tive effects. For a system of N elements with two states, moving along a periodic 
structure, we can introduce the distribution function p(a;, CTi, tT^v) for find- 
ing the particles i = 1..N in states ai — 1,2 with particle i at position x + si 
along a linear coordinate. This system thus becomes an effective 2^-state system 
described by 2^ equations 

JV 

dtp{x,ai,..,aN) + dJ{x,ai,..,aN) = - uj^j {x + is)p{x, ai,.., gn) 

i=l 

JV 

+ ^LJ^^{x + is)p{x,ai, ..,ai, ..,aN) (32) 

i=l 

Here, u}a{x) = Ucr{x+l) are the individual transition rates defined in the previous 
sections and the bar denotes the opposite state, i.e, 1 = 2,2 = 1. The currents 
are given by 

j{x, CTi, .., ctat) = [-ksTdxP - Npdxw{x, ai, .., vn) - Npfe^t] , (33) 
with the potential 

JV 

N 



w{x,ai,..,aN) = ^^Wai{x + is) , (34) 
=1 

defined as a sum of individual particle potentials. Here, rjN denotes the total 
friction which is assumed to scale linearly with N and /ext is the externally 
applied force per motor. 

In order to reduce the number of equations and to obtain a tractable descrip- 
tion, we introduce the average density of particles found at position ^ = a; mod /, 
< < I relative to the potential period: 

Pk{0=<Pk{0> , (35) 

where 

1 ^ 

Pk{0 = j;^J2^k,aAx + is-0 ■ (36) 

i=l 

Here, we have introduced the notation 

< a >= lim — / dx a{x,ai, ..,aN)p{x,ai, ..,aN) , (37) 
m-oo 2m J_^i ^f^^ 



for averages over the distribution p which we assume to be periodic, p{x, fJi, .., ctjv) = 
p{x + l,(Ti, ..jCtat) normalized over one period I, < 1 >= 1. The densities Pfe(^) 
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satisfy the normalization condition Eq. (20) and behave Hke a single particle 
two-state model 

dtPi + d^Ji = -LUiPi + LU2P2 

dtP2 + 5j J2 = wiPi - W2P2 , (38) 
however with the currents 

JkiO = -^d^Pk- < PkiOv > , (39) 

where v{x, txi, .., ctat) = —[dxw{x, txi, .., crjv) + /cxt]/??- From now on, we consider 
a large number A'' of motors and we assume that the period of motors is incom- 
mensurate with the potential period, s/l irrational. In this case, the total particle 
distribution function P1-I-P2 =< J2iLi S{^ — x — is)/N > becomes homogeneous: 

Pii0+P2{0 = j + Oil/N) , (40) 

and we can approximate < /3fe(0^ >=< pk{0 >< ^ > +0{l/N). Ignoring terms 
of order 1/A'' including the diffusive term in Eq. (39) the currents simplify to 



where 



=< V >= 



[ d^{Pld^Wi + P2d^W2) + /ext 
^0 



(41) 



(42) 



We have found a simple mean-field theory which is very useful to explore the 
properties of a many- motor system. Ignoring all corrections in 1/N, we finally 
obtain 



dtPi + vd^Pi = -(wi + (j2)Pi + 



L02 

T 



V 



i d^Pid^Aw + u, 

Jo 



(43) 
(44) 



which describes the time-evolution of Pi(^) = 1/Z — P2(C)- 



4.2 Steady states 

First, we look at the properties of steady state solutions with dfPi = and 
constant velocity which obey [46] 



vdfPi 



-{UJI + UJ2)Pl + ^2/1 



(45) 



where ii is a parameter. Eq. (45) can be solved analytically for simple choices of 
the transition rates. A more general approach is a power expansion of the steady 
state in the velocity 



(46) 



n=0 
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where the P^"^ satisfy the recursion relation P^"^ = -(9^Pi^"-^V(wi + ^2) 
with Pi^^ ~ ll)2/{uji + u>2)l- Inserting this solution in Eq. (44), we obtain a 
force- velocity relation /ext = f{v), with 

00 

/(z;) = 4°) + (^ + ^a))^ + ^^H^n ^ (47) 

n=2 

where ^ 

Fj^^ = - [ d^P^"-^ AW , (48) 







and the prime denotes a derivative with respect to ^. Using the definition Eq. 
(27) for the transition rates, we find the spontaneous force 



F^;^> = -td^ ^ (49) 







This force is zero for a system which is not chemically driven, i.e. i7 = or if 
the system is symmetric. Similarly, 



.(1) _ /' {AWre-^^/'-^-n'AW 
" ~ Jo Za;2(/2 + l + e-^WfeBT)3 ^^"^ 



is an effective friction coefficient which results from chemical transitions. For a 

passive system, this friction must be positive, F^2=o ^' "^hich is indeed the 
case. For =/= 0, however, /^^^ can become negative if fi' AW is positive. This 
criterion implies that the function 0{^) has maxima and minima at the same 
positions as AW . For Wi periodic and W2 constant this suggests e.g. placing 
the maxima of fi at the positions of minimal Wi which is the idea of localized 
excitations. 

The fact that f//^ can become negative suggests the possibility of an insta- 
bility in the system. This can be discussed most easily by first considering a 
symmetric system with symmetric potentials and transition rates. In this case, 
the system is invariant under (w, /cxt) (— — /ext), indicating that all even 
coefficients Z^^"'' = vanish by symmetry. In this case. 



f{v) = {n + F^^^)v + F^^\^ + 0{v') . (51) 



If Fq'' is negative and decreases with increasing amplitude of Q, an instability 
occurs as soon as Q = Qc for which 

FS = -n ■ (52) 

For Q > Qc, the function f(v) has a maximum, a minimum and a velocity 
interval where the effective friction 



(53) 
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becomes negative, see Fig. 8. If no external force is applied, the system does not 
move for Q < but two moving solutions bifurcate for Q > Qc with opposite 
velocity while the non-moving solution becomes unstable as revealed by a linear 
stability analysis. The system therefore a,t fi = [2^ undergoes a bifurcation with 
spontaneous symmetry breaking which allows for the existence of spontaneous 
motion in the symmetric case. 




Fig. 8. Relation between velocity v and externally applied force /oxt for a symmetric 
system. For i7 < Qc the system is passive. A dynamic transition with spontaneous 
symmetry-breaking occurs for Q = flc- For Q > Qc, the system is actively moving 
with velocity v± . The hysteresis is indicated by a broken line. 

If an external force is applied, the system can work against this force until 
the instability is reached at the minimum of the force velocity relation. Further 
increase of the force leads to a discontinuous change of the velocity. The system 
therefore shows the signature of a first-order transition with hystcrctic behavior 
and for /ext = bistability where the selection of one of the stable states depends 
on the history of the system. 

These arguments can be extended to the general case of a system with polar 
symmetry. In this case, all coefficients and in particular the spontaneous force 
Fq^ are nonzero. Dynamical transitions and instabilities do still exist but now 
occur typically for nonzero velocity and in an asymmetric way. A transition oc- 
curs if for increasing Q a point with ?7off = exists for a critical velocity Vc 
and for a critical force fc- For Q > He, a finite region with r/eff = emerges 
and the function f{v) has a maximum and a miniuimu. In this case, discontin- 
uous changes of the velocity occur as a function of applied force as soon as the 
force reaches the extremal values where the system becomes unstable and the 
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transition shows a hysteresis. A behavior with the characteristic as predicted by 
this theoretical analysis has been observed in motility assay experiments where 
electric fields are used to exert an external force [18]. 



4.3 Spontaneous oscillations 

Another important consequence of the collectivity is the possibility to generate 
oscillatory motion if the system is acting together with elastic elements, see Fig. 
7 (c). This corresponds to adding an elastic element to the force balance Eq. 
(44), which leads to 



1 

V = — 
V 



[ diP^AW + f,^,-kx 
Jo 



(54) 



where K = kN is the modulus of an external elastic element, k, the modulus 
per motor and 

x{t) = I dt'v{t') , (55) 
Jo 

is the total displacement of the system. First, we assume that the modulus 
k is small. In this case, the force varies slowly and we can use an adiabatic 
approximation assuming that at any given time the system is in one of the 
steady states with 

/ext -kx^ f{v) , (56) 

where f{v) is the force- velocity relationship introduced in Eq. (47). The change 
in velocity can be expressed as 

^ = -^ • (57) 

This equation demonstrates that we expect an oscillatory instability to occur in 
this adiabatic limit, exactly if 

77 + F«<0 . (58) 

In this case, the non-moving state with = is unstable. As soon as the system 
starts moving an elastic force is created which opposes motion and leads to a de- 
crease in velocity according to Eq. (57). As soon as an instability is reached with 
T]eS = 0, the velocity changes discontinuously in the opposite direction and the 
process is repeated. This scenario leads to characteristic relaxation oscillations 
of the position with a sawtooth like shape. For a symmetric system, the forward 
and backward parts of this process are identical and oscillations are symmetric 
with respect to f — > —t. In the asymmetric case, however, the asymmetry of the 
system is reflected in the shape of the oscillations, see inset of Fig. 9. 

This adiabatic approximation is valid as long as the characteristic time rjce/k 
is much longer than the transition times, i.e. rjes/k » ^ -l-w^^. For large k and 
in particular close to the discontinuous transition with T/eff = 0, the full dynamic 
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max 




no 

oscillations 



Fig. 9. Phase diagram for spontaneous oscillations as a function of the excitation ac- 
tivity i7 and the modulus per motor k of an external elastic element. The resting state 
becomes unstable with respect to oscillations along a line with O = Qc{k). The insets 
show typical examples for oscillatory motion (displacement versus time); onset of si- 
nusoidal motion at the instability and nonlinear oscillations with cusp-like shape for 
small k and large fi. 



equations have to be taken into account. The onset of oscillatory motion can be 
determined by a linear stability analysis of the non-moving state with w = 0, 
and X = Xq = —f'o'/k. For a small deviation from this initial state we write the 
ansatz 



v{t) = ue** , 



(59) 
(60) 



where R = P^"^ and s = iui + t is complex. Here, u) is an oscillation frequency 
and r a relaxation time. Inserting these expressions in Eq. (43), we find to linear 
order 

where we have defined a{x) = ui{x) + co2{x). Possible values of s are solutions 
to 

k f' , R'AW 
V+- = - d^^—- , (62) 

s Jo a + s 

As long as T < 0, the initial state is locally stable, it becomes unstable for r = 
towards oscillatory motion with frequency co. This instability is determined by 



Force and Motion Generation of Molecular Motors: A Generic Description 



23 



setting s = iujc- The frequency Qc at the instability obeys 

^ = - f di^^^BlAW (63) 

rl ^2 

k= / d^^-^R'AW . (64) 

Jo ^2+^2 

For fc = 0, we recover the criterion (58) with = 0. The frequency at the 
instability increases as Uc ~ k^^'^ for increasing values of k. However, for k > fcmax 
with 

kmax < - I RAW" , (65) 



^0 

no instability occurs for arbitrarily large SI see Fig. 9. Note, that fcmax is of the 
order of magnitude of the average negative curvature of the potential AW and 
can be roughly estimated by fcmax — U/l"^, where U is the potential period. Note, 
that since k a modulus per motor, a motor collection of N motors can therefore 
undergo oscillatory motion if working against an elastic element with rigidity 
smaller than K < kmax^- 

The nonlinear relation between force and velocity of steady states can be 
generalized to time-dependent external forces. Assuming that a limit cycle exists 
which is characterized by the Fourier representations 

oo 

/ext(0= E /"e"* 
n= — oo 
oo 

v{t) = E ^"^ 



iujt 



(66) 

a generalization of the expansion (47) is given by 

f„ = Fi°^+Fi]jvk+Fi'lviVrn + Fi%^VlVrnVo + Oiv') . (67) 

The spontaneous force Fn^"' = F^'^^Sno is nonvanishing only for the time inde- 
pendent mode n = 0. The linear response coefficient 



_k_ 
^ iwl 



,R'AW' 
I d£ ^ 

Jo Q! -I- ZLOl 



(68) 



vanishes at the instability. The real and the imaginary part of Fi^ii corresponds 
to an effective friction r]es and to an effective elastic modulus kes/ico, respec- 
tively. Because of the non equilibrium nature of the systems both of these effec- 
tive linear response coefficients can become negative which is impossible in an 
equilibrium situation. 
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Higher order coefficients can be obtained in a systematic way [47]. The ex- 
pansion (67) characterizes the history-dependent response of the system to an 
external force which is periodic in time for small velocity amplitudes. As a re- 

(n) 

suit of time translation symmetry, the coefficients obey -F„j ; ^ Sm,ii+..+i„- 
Spontaneous oscillations are solutions to Eq. (67) with /„ = for all n which 
are well described near the instability. Far away from the instability nonlinear- 
ities become dominant. In the regime of small k, they are well captured by the 
adiabatic approximation given by Eq. (57). 

5 Discussion and outlook 

In the previous sections, we have described a generic framework to model the 
force and motion generation of molecular motors. The coupling of a chemical 
reaction to internal conformations and spatial degrees of freedom of a motor 
enzyme is described by an overdamped stochastic dynamics. Reducing the de- 
scription to a simple two state model has the advantage to be sufficiently simple 
for a theoretical analysis but still flexible to allow for a description of conditions 
of operation similar to those of motor molecules such as myosin or kinesin. 

A quantitative modelization of motor enzymes on the molecular level is diffi- 
cult. One the one hand experimental data is limited, on the other hand simplified 
models are not designed to describe all details of the functioning of complex en- 
zymes. The properties of the two-state model discussed here depend of the choice 
of the potential shapes and on the conformation-dependence of the transition 
rates described by the functions a{x) and P{x). These functions cannot be de- 
duced from the molecular structure and are difficult to obtain experimentally. 
This problem is often referred to as the unknown strain-dependence of transi- 
tions [37, 43]. The known force- velocity relation for single kinesin molecules is 
close to linear [41, 17]. Different models which have been used to describe kinesin 
motion, all produce almost linear force-velocity relations and describe well the 
observed orders of magnitude [26, 38, 43, 44] which demonstrates that the basic 
properties of force-generation do not depend much on the details of a chosen 
model. For these reasons, we refrain from a full quantitative comparison and 
focus our discussion on physical principles. 

In biological situations, molecular motors are involved in processes such as 
muscle contraction, cell division and flagellar beating which are complex and 
involve large numbers of motors together with other enzymc!s which control and 
regulate the function of these systems. A physical approach to more complex 
situations is to avoid the complexity of control systems and to investigate the 
types of behavior which can result from the activity of many motors alone. 
The properties of larger-scale active systems can be studied using concepts of 
statistical physics. An important advantage of two-state models in this context is 
the possibility to start from a simple model for individual motors and to address 
the properties of many-motor systems. From such an approach the possibility 
of instabilities, dynamic transitions and spontaneous oscillations follows quite 
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naturally. We therefore predict such behaviors to occur in biological situations. 
Related collective phenomena have recently been shown to exist also in other 
models [48, 49] which shows that the existence of dynamic instabilities in this 
class of systems is quite general. The detailed conditions for which dynamic 
transitions will occur in biological systems such as coupled myosins or kinesins 
should depend on various parameters such as salt concentrations, temperature, 
motor density etc. Our approach docs not predict for which situations these new 
behaviors can be found but provides a classification of the possible self-organized 
behaviors of such systems. 

In vitro motility assays which use purified motors and filaments in an artificial 
environment are ideal systems to test these ideas and to study motor action in the 
absence of further biochemical regulatory systems. Motors grafted to a substrate 
set fluorescently marked filaments in motion which can easily bo observed by 
standard light microscopy. The force-velocity relationship of the actin-myosin 
system has recently been determined by constraining filament motion along one 
dimension by using linear grooves of 1/im in diameter as a substrate [18, 50], see 
Fig. 3 (d). Electric fields E applied parallel to the aligned and negatively charged 
actin filaments induce external forces which are homogeneously distributed along 
the filaments. The electrophoretic mobility /i = v/E, where v is the velocity 
induced by an electric field E, measured for free filaments can be used to obtain 
an estimate for the external force per unit length / ~ 2tt iirjw / \n{r / d) E acting 
eflfectively on the filament. Here, eta^ is the viscosity of the solution and r and 
d are the radius of the filament and the distance between filament axis and the 
substrate. With this argument, we find that a field of E c^^ lO^V/m induces a 
force per length of / ~ IpN/iim which is the order of magnitude of stall forces 
of an actin-myosin system [18] . 

In addition to providing the first force- velocity relations for a purified actin- 
myosin system, these experiments have revealed unexpected behavior for high 
myosin densities close to stall conditions. With increasing force applied against 
the direction of motion, filaments slow down, however before stalling completely 
an abrupt change in velocity to reverse motion is observed. Since a histogram 
of the velocity distribution reveals two distinct maxima for diS^erent coexist- 
ing velocities and discontinuous changes of the velocity occur with a hysteretic 
behavior, a natural interpretation of this observation is a dynamic first order 
transition of the type discussed in the previous sections. However, experiments 
of this type are difficult and measured values have significant uncertainties; more 
data is necessary to clearly prove the existence of the transition. Interestingly, 
macroscopic measurements of the behavior of muscles show near stalling con- 
ditions an unexpected rapid drop in velocity and a contracting state which is 
difficult to stabilize experimentally by feedback systems [51]. This behavior could 
result from the averaged mechanical properties of many different filaments which 
individually show a transition of the type observed in motility assays. 

In recent years, an increasing number of situation have been observed for 
which the concepts of dynamic instabilities in many-motor systems could be 
relevant. It is well known that the flight muscles of many insects generate os- 
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cillations themselves. These so-called asynchronous muscles show oscillatory be- 
havior which is not triggered by a periodic nerve-signal [52]. These systems 
are complete muscles for which the role of biochemical control systems could 
be important. Recently, it has become possible to study the active mechanical 
properties of single myofibrils, i.e. contractile fibers within muscle cells, by using 
microneedles. Ordinary muscle myofibrils which under normal conditions simply 
contract, show spontaneous oscillations if put in specific conditions such as in- 
creased ADP concentration. Both tension- and length- oscillations were observed 
[53]. It was demonstrated that such oscillations continue to exist if regulatory 
enzymes bound to the actin- myosin motors were removed [54]. These observa- 
tions demonstrate the general possibility of oscillations in many muscles. The 
facts that in muscular structures, which couple contractile elements elastically, 
oscillations arc typically observed while the motility assay which lacks elastic 
elements reveals a first-order transition is fully consistent with the collective 
effects discussed here. 

Another important type of oscillating many-motor systems arc cilia and 
flagella. These are hair-like appendages of many cells which are used for self- 
propulsion and to stir the surrounding fluid. They share the characteristic archi- 
tecture of their core structure, the axoneme, a common structural motive that 
was developed early in evolution. It is characterized by nine parallel pairs of mi- 
crotubules which are arranged in a circular fashion together with a large number 
of dynein molecular motors. In the presence of ATP, the dyncins attached to the 
microtubules generate relative forces while acting on neighboring microtubules. 
The resulting internal stresses induce the bending and wave-like motion of the 
axoncmc [1]. A simple two-dimensional model of filaments driven internally by 
molecular motors can be used to demonstrate that a dynamic instability of the 
many-motor system most naturally generates oscillating motion and wave-like 
propagating shapes which can be used for self-propulsion [55]. This idea that 
simple patterns of motion are induced via a dynamic instability is supported by 
the observations that flagellar dyneins are able to generate oscillatory motion 
on microtubules [56], and that isolated and de-membranatcd flagella in solution 
containing ATP above a threshold concentration swim with a simple wave-like 
motion [57]. This suggests, that basic wave-like patterns of flagellar motion are 
generated by a self-organized motor-induced mechanism and that sophisticated 
regulatory systems could have evolved at a later time to fine tune the system 
and to generate new patterns of motion. 

The author thanks A. Ajdari, M. Badoual, M. Bornens, S. Camalct, P. 
Chaikin, R. Everaers, A. Maggs, A. Ott, A. Parmeggiani, J. Prost, D. Riveline 
and C. Wiggins for stimulating interactions. 
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